３．研究系及び研究施設の現状
3-1 論文発表状況
3-1-1 論文の発表状況
分子研では毎年 A nnual R eview（英文）を発刊し，これに発表した全ての学術論文のリストを記載している。
論文の発表状況
間 期 象 対 集
編 ANNUALREVIEW 原著論文の数 総説等の数
∼1978.8. 1978 25 13 .
9 . 8 7 9
1 ∼1979.8. 1979 55 7 .
9 . 9 7 9
1 ∼1980.8. 1980 85 21 .
9 . 0 8 9
1 ∼1981.8. 1981 114 24 .
9 . 1 8 9
1 ∼1982.8. 1982 149 14 .
9 . 2 8 9
1 ∼1983.8. 1983 177 29 .
9 . 3 8 9
1 ∼1984.8. 1984 153 26 .
9 . 4 8 9
1 ∼1985.8. 1985 196 31 .
9 . 5 8 9
1 ∼1986.8. 1986 207 45 .
9 . 6 8 9
1 ∼1987.8. 1987 287 42 .
9 . 7 8 9
1 ∼1988.8. 1988 247 39 .
9 . 8 8 9
1 ∼1989.8. 1989 281 60 .
9 . 9 8 9
1 ∼1990.8. 1990 320 60 .
9 . 0 9 9
1 ∼1991.8. 1991 260 23 .
9 . 1 9 9
1 ∼1992.8. 1992 303 41 .
9 . 2 9 9
1 ∼1993.8. 1993 298 41 .
9 . 3 9 9
1 ∼1994.8. 1994 211 26 .
9 . 4 9 9
1 ∼1995.8. 1995 293 23 .
9 . 5 9 9
1 ∼1996.8. 1996 332 40 .
9 . 6 9 9
1 ∼1997.8. 1997 403 41 .
9 . 7 9 9
1 ∼1998.8. 1998 402 44 .
9 . 8 9 9
1 ∼1999.8. 1999 401 47 .
9 . 9 9 9
1 ∼2000.8. 2000 337 30 .
9 . 0 0 0
2 ∼2001.8. 2001 405 65 .
9 . 1 0 0
2 ∼2002.8. 2002 489 59
3-1-2 論文の引用状況
論文の引用状況については，一昨年度
1)
及び昨年度
2)
の「分子研リポート」で報告したが，本年度も続けて報告することにする。 論文の引用数については，米国 IS I 社（T he Institute for S cientific Information）
3)
の引用統計データベースに基づく調査が標 準になりつつある。しかし，このデータベースを利用しても，新聞や雑誌などで公開されている多くの研究機関のランキング結果を 見るとき，総論文数や総被引用数のように，研究機関の研究者の数に強く依存する量を基準としている場合がほとんであり，それ らの数値は研究者が多いということを意味しているに過ぎず，研究レベルの本質を見ているとは言えない。例えば，IS I社の日本
支社が発表したホームページの資料
4)
では，総引用数を元に日本のトップ１０機関をランク付けしているが，上位に大規模な大学 である，所謂旧帝国大学等が名前を連ね，ランキングの結果と研究者の数に強い相関があることは一目瞭然である。この問題点 については，以前の分子研リポート
1),2)
においても議論した。例外は，国立情報学研究所の根岸正光教授らの考察
5)
で，そこでは， 論文１報あたりの平均被引用数（引用度）を基準とすべきであるとしている。この量は研究者の数に依存しないので，「研究の質」 についての大変信頼性の高い指標を与える。本リポートでも，引用度を指標とすることにする。
文献４と５が今年度発表された最新の結果であるので，それらについて，まず，まとめることにする。すなわち，昨年の分子研リポー トの結果
2)
をアップデートする。文献５では，ISI社のデータベース中の論文の所属機関のうち，大学・大学共同利用機関等（「大学 等」と呼ぶことにする）に所属するものを調査の対象としている。以下，本リポートでも大学等に対象を絞ることにする。上でも述べ たように，文献４はIS I社の日本支社が発表したホームページの資料である。そこでは，１９９２年１月から２００２年１２月までの１１年間 に発表された論文のISI社のデータベース中の論文の総引用数を元に日本のトップ１０機関をランク付けしている。このホームペー ジの表では，総引用数のほかに，総論文数と引用度（論文１報あたりの平均被引用数）も掲載されているので，引用度によるラン ク付けを表１にまとめた（分子研がランクインしている化学分野における結果である）。昨年の結果から調査期間が一年ずれただ けなので，４位と５位の名古屋大学と京都大学の順位が入れ替わっただけで，他には変化がなかった。分子研に関しては相変わ
表１日本の大学等の分野別論文引用度 表２ 日本の大学等の分野別論文引用度指数 分野：化学（１９９２―２００２） 分野：化学（１９９２―２００１）
位
順 大 学 等 論文引用度 順 位 大 学 等 引用度指数 1 分子科学研究所 11.33 1 岡崎国立共同研究機構 165 2 東京大学 9.38 2 東京大学 135 3 北海道大学 9.07 3 名古屋大学 126 4 名古屋大学 8.98 4 京都大学 125 5 京都大学 8.71 5 北海道大学 124 6 東京理科大学 8.59 6 東京都立大学 121 7 九州大学 8.33 7 東京理科大学 118 8 大阪大学 8.22 8 大阪大学 115 9 東京工業大学 7.99 9 九州大学 113 0
1 東北大学 7.71 10 東京工業大学 110 11 東北大学 108
2
1 早稲田大学 107 3
1 千葉大学 105 4
1 大阪市立大学 102 5
1 広島大学 100
らず第１位の地位を守っていることが判明した。特に，分子研が２位に1.95ポイントの差をつけているの対し，２位から１０位までの 間に1.67ポイントの差しかないことに注目されたい。文献４の資料では，まず，総引用数がトップ１０の研究機関までで「足切り」して いるので，研究者の数が多い大学に比べて少ない分子研は他の分野（例えば，物理）では引用度が上位にありながら，考慮から はずれている（それにもかかわらず，化学の分野では引用度の圧倒的な高さによって研究者の数の少なさをカバーして，ランクイ ンしたということもできる）。
今年度得られた新しいデータの二つ目は根岸氏による文献５である。今回は１９９２年から２００１年までの１０年間についての解析 である。ここでは，引用度指数という量でランク付けがされている。引用度指数については，前回詳しく説明した
2)
が，ここでは，引 用度指数が１００の研究機関は，その分野で我が国で平均的な引用度の論文を出しており，２００ならば平均の２倍の引用度になっ ていることに注意すれば十分であろう。引用度指数による化学における分野別ランク付けを表２にまとめた。
5)
考慮する研究機関 は，論文数上位１５機関としている。
5)
前年度の同様の考察では，引用数合計の上位５０位以内の研究機関（但し，論文数が１０件以 上）を対象としたので，足切りの規準が変わったことに注意されたい（そのために，前回化学で２位にランクインした総研大が，今 回はランク外となった）。化学では分子研が圧倒的に全国第１位であることが再び確認できる。
以上，分子科学研究所が（特に化学の分野で）我が国で極めて高い位置づけを得ていることを示した。最後に，国際的な比較 を試みることにする。これについて，文献６と７が見つかったので，以下にその内容について説明する。
文献６では，理工系分野と生物・医学系分野について，それぞれ，引用度を日本の大学等についてまとめた（表３と表４）。このよ うに大きく分野を２つに分けて考慮するのは，これら二大研究分野の習慣によって引用の仕方に相当の差があるためである。例え ば，分子生物学・遺伝学の分野の平均引用度は17.1なのに対し，工学の分野では3.0である。
5)
文献６によれば，理工系分野では， 岡崎国立共同研究機構の分子科学研究所が我が国第１位，生物・医学系では，岡崎国立共同研究機構の基礎生物学研究所と 生理学研究所がそれぞれ我が国第１位と第２位という結果となった。これらの表に，参考として，世界の一流大学の数値が示され
表３日本の大学等の論文引用度 表４ 日本の大学等の論文引用度 分野：理工系（１９８１―１９９１） 分野：生物・医学系（１９８１―１９９１） 位
順 研究機関 論文引用度 順 位 研究機関 論文引用度 1 分子科学研究所 9.89 1 基礎生物学研究所 20.74 2 東京大学 8.22 2 生理学研究所 16.74 3 京都大学 7.05 3 神戸大学 13.80 4 大阪大学 6.54 4 宮崎医科大学 12.77 5 東京工業大学 6.49 5 国立遺伝学研究所 11.83 6 名古屋大学 6.29 6 静岡薬科大学 11.18 7 東北大学 6.27 7 大阪大学 10.02
8 自治医科大学 9.97
考
参 高等研究所(プリンストン) 17.47 9 京都大学 9.96 ハーバード大学 15.71 11 東京大学 9.87 12 東京工業大学 9.39 13 北海道大学 8.86
考
参 カリフォルニア工科大学 24.49
表５ 世界の大学等の分野別論文引用度 分野：化学（１９９１―２０００）
位
順 大 学 等 論文引用度 1 ハーバード大学 21.6 2 カリフォルニア工科大学 20.9 3 スタンフォード大学 17.9 4 マサチューセッツ工科大学 16.5 5
ア ニ ル ォ フ リ カ
， 学 大 ル ネ ー コ
校 ー レ ク ー バ 学 大
4 . 5 1 7 ミネソタ大学 14.9 8 イリノイ大学 13.3 9 テキサス大学 13.0 0
1 スイス工科大学チューリッヒ校 11.9 1
1 ケンブリッジ大学 11.7 2
1 ウィスコンシン大学 11.4 3
1 パデュー大学 11.3 4
1 ミュンヘン工科大学 10.9
ているが，理工系では分子科学研究所が世界の一流に対して数ポイントの差をつけられているものの，大いに健闘していること が分かる。また，生物・医学系では，基礎生物学研究所は世界の一流と比べて遜色ないと言え，生理学研究所は数ポイントの差 以内で健闘していることが分かる。
文献７は更に新しいデータであり，ケンブリッジ大学の化学科が出版している雑誌の資料である。結果を表５にまとめた。そこで は，総被引用数の世界ランキング２０位までを対象としているので，研究者の多くない分子研はランクからはずれている。しかし，表 ５と表１を比べると，化学の分野で我が国第１位の分子研が世界では大体１１位∼１３位ぐらいで，ケンブリッジ大学やウィスコンシン
大学とほぼ同レベルであることが分かった。
最後にまとめると，引用度は，（常に顔を合わして交流を持っている）有力研究者の数が多く，また，多くの有力雑誌が刊行され ている欧米（特に米国）における研究機関が圧倒的に有利な立場にあることを考慮すると，岡崎国立共同研究機構の三研究所 の世界順位は非常に高いものであると言える。
（分子基礎理論第一研究部門 岡本祐幸 記）
参考文献
1. 分子研リポート2001, pp. 62–66. 2. 分子研リポート2002, pp. 74–77. 3. http://www.isinet.com/
4. http://www.isinet.com/japan/news/20030408.html
5. 根岸正光 , 「大学ランキング 2004」, 朝日新聞社 , pp. 156–163 (2003). 6. A. Anderson et al., “Science in Japan,” Science 258, 561 (1992). 7. C hem@ C am (Issue 13, W inter 2002, C ambridge University), p. 7.
3-2 理論研究系
分子基礎理論第一研究部門
永 瀬 茂（教授）
A -1)専門領域：理論化学、計算化学
A -2)研究課題：
a) 高周期元素の特性を利用した分子の設計と反応 b)分子の立体的な形と大きさを利用した分子設計と反応 c) ナノスケールでの分子設計理論と計算システム
A -3)研究活動の概略と主な成果
a) 周期表には利用できる元素は８０種類以上もあり，これらの組み合わせから生まれる新しい結合と相互作用は機能発 現に無限の可能性を秘めている。高周期元素を骨格にもつ新規な多重結合分子，芳香族分子，クラスター等の開拓と 構築を行い，特異な構造と電子状態および機能を実験と密に連携して理論的に明らかにした。
b)分子の特性は，元素の組み合わせばかりでなく，立体的な形状とサイズおよび柔軟さにも大きく支配される。サイズ の大きな分子が作る外部空間および内部空間は新しい機能発現のための相互作用場として利用できる。このために， フラーレン，ナノチューブ，カーボンピーポッド，ホストーゲスト超分子，酵素反応などの系統的な理論研究をおこ なった。また分子デバイスへの応用のために，球状分子に内部にとりこめられた原子や分子の動的挙動の制御を化 学修飾により行った。
c) 実験で興味ある系を広く対象とするには，分子理論と計算の適用範囲を大きく拡張する必要がある。このとき一番 の問題は，取り扱う系の大きさが大きくなると計算時間が急激に増加することである。この問題を解決するために 有用な分子理論と計算法および高速並列システムの開発を行っている。
B -1) 学術論文
K. KOBAYASHI and S. NAGASE, “Theoretical Calculations of Vibrational Modes in Endohedral Metallofullerenes: La@C82 and Sc2@C84,” Mol. Phys. (a special issue for Prof. Yoshimine) 101, 249–254 (2003).
G. M. A. RAHMAN, Y. MAEDA, T. WAKAHARA, M. KAKO, S. SATO, M. OKAMURA, T. AKASAKA, K. KOBAYASHI and S. NAGASE, “Photochemical Bissilylation of C70 with Disilane,” ITE Lett. Batt. New Tech. Med. 4, 60– 66 (2003).
T. WAKAHARA, G. M. A. RAHMAN, Y. MAEDA, M. KAKO, S. SATO, M. OKAMURA, T. AKASAKA, K. KOBAYASHI and S. NAGASE, “Redox Properties of Carbosilylated and Hydrosilykated Fullerene Derivatives,” ITE Lett. Batt. New Tech. Med. 4, 67–73 (2003).
T. WAKAHARA, M. KAKO, Y. MAEDA, T. AKASAKA, K. KOBAYASHI and S. NAGASE, “Synthesis and Characterization of Cyclic Silicon Compounds of Fullerenes,” Cur. Org. Chem. 7, 927–943 (2003).
K. NAGAYOSHI, K. KITAURA, S. KOSEKI, S. RE, K. KOBAYASHI, Y. -K. CHOE and S. NAGASE, “Calculation of Packing Structure of Methanol Solid Using Ab Initio Lattice Energy at the MP2 Level,” Chem. Phys. Lett. 369, 597–604 (2003).
J. LU, S. RE, Y. -K. CHOE, S. NAGASE, Y. ZHOU, R. HAN, L. PENG, X. ZHANG and X. ZHAO, “Theoretical Identification of Carbon Clusters C20: Prevalence of the Monocyclic Isomer and Existences of the Smallest Fullerene and Bowl Isomer,” Phys. Rev. B 67, 125415 (7 pages) (2003).
J. LU and S. NAGASE, “Structural and Electronic Properties of Metal-Encapsulated Silicon Clusters in a Large Size Range,” Phys. Rev. Lett. 90, 115506 (4 pages) (2003).
J. LU and S. NAGASE, “Metal-Doped Germanium Clusters MGen at the Sizes of n = 12 and 10: Divergence of Growth Patterns from the MSin Clusters,” Chem. Phys. Lett. 372, 394–398 (2003).
Z. SLANINA, K. KOBAYASHI and S. NAGASE, “Ca@C72 IPR and Non-IPR Structures; Computed Temperature Development of their Relative Concentrations,” Chem. Phys. Lett. 372, 810–814 (2003).
N. TAKAGI, K. YAMAZAKI and S. NAGASE, “Theoretical Investigation of Triple Bonding between Transition Metal and Main Group Elements in (η5-C5H5)(CO)2M≡ER (M = Cr, Mo, W; E = Si, Ge, Sn, Pb; R = Terphenyl Groups),” Bull. Korean Chem. Soc. (a special issue) 24, 832–836 (2003).
K. NAGAYOSHI, T. IKEDA, K. KITAURA and S. NAGASE, “Computational Procedure of Lattice Energy Using the Ab Initio MO Method,” J. Theor. Comput. Chem. 2, 233–244 (2003).
K. KOBAYASHI, S. NAGASE, Y. MAEDA, T. WAKAHARA and T. AKASAKA, “La2@C80: Is the Circular Motion of Two La Atoms Controllable by Exohedral Addition? ” Chem. Phys. Lett. 374, 562–566 (2003).
K. KOBAYASHI, S. NAGASE and K. -P. DINSE, “A Theoretical Study of Spin Density Distributions and Isotropic Hyperfine Couplings of N and P atoms in N@C60, P@C60, N@C70, N@C60(CH2)6, and N@C60(SiH2)6,” Chem. Phys. Lett. 377, 93–98 (2003).
Y. MAEDA, G. M. A. RAHMAN, T. WAKAHARA, M. KAKO, M. OKAMURA, S. SATO, T. AKASAKA, K. KOBAYASHI and S. NAGASE, “Synthesis and Characterization of Tetrakis-silylated C60 Isomers,” J. Org. Chem. 68, 6791– 6794 (2003).
T. WAKAHARA, Y. MAEDA, M. KAKO, T. AKASAKA, K. KOBAYASHI and S. NAGASE, “Silylation of Fullerenes with Active Species in Photolysis of Polysilane,” J. Organomet. Chem. 685, 177–188 (2003).
T. TAJIMA, K. HATANO, T. SASAKI, T. SASAMORI, N. TAKEDA, N. TOKITOH, N. TAKAGI and S. NAGASE,
“Syntheses and Structures of Silicon Analogues of Cyclopropabenzenes,” J. Organomet. Chem. 686, 118–126 (2003). M. T. H. LIU, Y. -K. CHOE, M. KIMURA, K. KOBAYASHI, S. NAGASE, T. WAKAHARA, Y. NIINO, M. ISHITUKA and T. AKASAKA, “The Effect of Substituents on the Thermal Decomposition of Diazirines: Experimental and Computational Studies,” J. Org. Chem. 68, 7471–7478 (2003).
J. LU, S. NAGASE, L. PENG and S. ZHANG, “Strongly Size-Dependent Electronic Properties in C60-Encapsulated Zigzag Nanotubes and Lower Size Limit of Carbon Nanopeapods,” Phys. Rev. B 68, 121402 (4 pages) (2003).
T. WAKAHARA, Y. MATSUNAGA, A. KATAYAMA, Y. MAEDA, M. KAKO, T. AKASAKA, M. OKAMURA, T. KATO, Y. -K. CHOE, K. KOBAYASHI, S. NAGASE, H. HUANG and M. ATA, “A Comparison of the Photochemical Reactivity between N@C60 and C60: Photolysis with Disilirane,” Chem. Commun. 2940–2941 (2003).
Z. SLANINA, K. KOBAYASHI and S. NAGASE, “Temperature Development in a Set of C60H36 Isomers,” Chem. Phys. Lett. 382, 211–215 (2003).
Z. CHEN, A. HIRSCH, S. NAGASE, W. THIEL and P. v. R. SCHLEYER, “Spherical Sila- and Germa-Homoaromaticity,” J. Am. Chem. Soc. 125, 15507–15511 (2003).
B -3) 総説、著書
T. WAKAHARA, Y. MAEDA, S. OKUBO, J. KOBAYASHI, M. KONDO, T. AKASAKA, K. KOBAYASHI , S. NAGASE, T. KATO, K. YAMAMOTO and K. M. KADISH, “Endohedral Fullerene Ions: Synthesis, Structure and Reaction,” in Structure and Electronic Properties of Molecular Nanostructures, H. Kuzmany, J. Fink, M. Mehring and S. Roth, Eds., 43–46 (2002).
Z. SLANINA, K. KOBAYASHI and S. NAGASE, “Computing Temperature Development of Five Isomers of Ca@C72,” in Fullerenes and Nanotubes: The Building Blocks of Next Generation Nanodevices, D. M. Guldi, P. V. Kamat, and F. D’Souza, Eds., The Electrochemical Society, Inc.; Pennington, NJ, 13, 582–590 (2003).
永瀬 茂 , 「典型元素化学の特徴」, 有機金属反応剤ハンドブック, 化学同人 , 4–10 (2003). 永瀬 茂 , 「理論と計算による分子設計」, 先端化学シリーズ IV , 丸善 , 16–20 (2003).
B -4) 招待講演
S. NAGASE, “Endohedral Metallofullerenes: Cage Structures and Metal Motion,” 10th Korea-Japan Joint Symposium on Theoretical/Computational Chemistry, Pohang (Korea), January 2003.
K. KOBAYAHI, “Theoretical Prediction for Endohedral Metallofullerenes,” Technische Universität Darmsadt, Germany, May 2003.
K. KOBAYAHI and S. NAGASE, “Is the Motion of La Atoms inside C80 Controllable by Exohedral Addition? ” The 203th International Meeting of Electrochemical Society, Parris (France), April 2003.
永瀬 茂 , 「ナノサイズ分子の計算化学」, 第２回筑波大学スーパーコンピュータワークショップ , 筑波 , 2003 年 5月 . 永瀬 茂 , 「分子理論と計算化学の進展」, 秋田地区講演会（日本化学会東北支部）, 秋田 , 2003 年 11月 . 永瀬 茂 , 「高周期元素の多重結合の進展」, 第３０回ヘテロ原子化学討論会 , 富山 , 2003年 12月 .
小林 郁 , 「化学修飾によるフラーレン内部の金属原子の状態変換」, 日本物理学会（シンポジウム）, 岡山 , 2003 年 9月 .
B -7) 学会および社会的活動 学協会役員、委員
W A T OC（世界理論化学学会）委員 . 学会の組織委員
「10th Korea-Japan Joint Symposium on Theoretical/Computational Chemistry」共催 . 分子研研究会「固体表面における非熱的電子励起状態の化学」主催 .
A PA C T C 実行委員 .
分子構造総合討論会運営委員会幹事 . フラーレン・ナノチューブ研究会幹事 .
フラーレン若手の会世話人代表（小林 郁）. 学会誌編集委員
Silicon Chemistry, Subject Editor.
B -8) 他大学での講義、客員
千葉大学理学部化学科 , 集中講義「計算機有機化学」, 2003年 8 月 . 筑波大学先端学際領域研究センター併任教授，2002年 11 月 - . 筑波大学 T A R A センター, 客員研究員 , 2002年 1月 - .
C ) 研究活動の課題と展望
新素材開発において，分子の特性をいかにしてナノスケールの機能として発現させるかは最近の課題である。このために， 炭素を中心とする第２周期元素ばかりでなく大きな可能性をもつ高周期元素およびナノ構造の特性を最大限に活用する分 子の設計と反応が重要である。サイズの大きい分子はさまざまな形状をとれるので，形状の違いにより電子，光，磁気特性ば かりでなく，空孔の内径を調節することによりゲスト分子との相互作用と取り込み様式も大きく変化させることができる。これら の骨格に異種原子や高周期元素を加えると，変化のバリエーションを飛躍的に増大させることができる。ナノスケールでの 分子設計理論とコンピューターシミュレーション法を確立し，高い分子認識能をもつナノ分子カプセル，機能性超分子，疑似 タンパク質，デンドリマーおよび伝導性共役高分子を開発する。これらの分子を効率的に合成実現するためには，従来のよ うに小さい分子から順次組み上げていくのではなく，自己集合的に一度に組織化する機構の解明と理論予測はきわめて重 要である。また，現在の量子化学的手法は，小さな分子の設計や構造，電子状態，反応を精度よく取り扱えるが，ナノスケー ルでの取り扱いには飛躍的な進展が望まれている。
岡 本 祐 幸（助教授）
A -1)専門領域：生物化学物理、計算科学
A -2)研究課題
a) 蛋白質分子の第一原理からの立体構造予測問題および折り畳み問題 b)生体分子以外の系への拡張アンサンブル法の適用
A -3)研究活動の概略と主な成果
a) 蛋白質は自然界に存在する最も複雑な分子である。よって，その立体構造を予測することは（その生化学的機能との 関係上，極めて重要であるにもかかわらず）至難の業である。特に，理論的に第一原理から（自由エネルギーを最小化 することにより）立体構造を予測することは不可能と広く信じられている。それは，溶媒の効果を取り入れるのが困 難であるばかりでなく，系にエネルギー関数の極小状態が無数に存在するために，シミュレーションがそれらに留 まってしまって，世界最速のスーパーコンピューターをもってしても，最小エネルギー状態に到達するのが絶望的 であるからである。我々はシミュレーションがエネルギー極小状態に留まらない強力な計算手法を，蛋白質の立体 構造予測問題に適用することを提唱してきた。具体的には，徐冷法（ si mul ated anneal i ng ）及び拡張アンサンブル法
（generalized-ensemble algorithm）を導入し，これらの手法が小ペプチド系において従来の方法よりはるかに有効であ ることを示してきた。拡張アンサンブル法では，非ボルツマン的な重み因子に基づいて，ポテンシャルエネルギー空 間上の酔歩を実現することによって，エネルギー極小状態に留まるのを避ける。この手法の最大の特徴は唯一回の シミュレーションの結果から，最小エネルギー状態ばかりでなく，物理量の任意の温度におけるアンサンブル平均 を 求 め る こ と が で き る こ とで あ る 。拡張 ア ン サ ンブ ル 法 の 代表 的 な 例 がマ ル チ カ ノニ カ ル 法（ mul ti c anoni c al algorithm）と焼き戻し法（simulated tempering）であるが，これらの二手法ではその重み因子を決定することが自明で はない。この問題を克服するため，我々は新たにT sallis統計に基づく拡張アンサンブル法やレプリカ交換法（replica- exchange method）の分子動力学法版を開発したりしてきた。特に，レプリカ交換分子動力学法はその適用が簡便であ るために，我々の発表とともにすぐに受け入れられ，現在では国内外のタンパク質の折りたたみシミュレーション における有力グループが相次いで採用している。更には，正確な溶媒の効果をエネルギー関数に取り入れていくこ とも大切であるが，距離に依存した誘電率で表すもの（レベル１）や溶質の溶媒への露出表面積に比例する項（レベ ル２）を試すとともに，厳密な溶媒効果(レベル３）として，R IS M や S PT などの液体の統計力学に基づくものや水分 子をあらわにシミュレーションに取り入れること等を検討してきた。
本年度は，２０００年に我々が開発したレプリカ交換マルチカノニカル法（R E MUC A ）とマルチカノニカルレプリカ交 換法（ MUC A R E M ）の有効性を，小ペプチド系のモンテカルロシミュレーションで確固としたものであることを示 した。そして，この２手法がタンパク質の折り畳みシミュレーションに特に適していることを示した（Mitsutake, Sugita and Okamoto, J. Chem. Phys. (2003)）。これらの新しい拡張アンサンブル法をレベル３の厳密な溶媒効果を取 り入れた（T IP3Pの水分子をあらわに取り入れた）アミノ酸数が十数個の小ペプチド系に適用することによって，広 く使われているA MB E R ，C HA R MM，OPL S，GR OMOSなどの標準的なエネルギー関数（力場）が蛋白質の立体構造予 測が可能な程の精度を持つか否かを調べてきたが，この判定に今年度やっと一つの区切りをつけることができた
（Yoda, Sugita and Okamoto, Chem. Phys. Lett. (2004)）。この判定には，エネルギー極小状態に留まらず，広く構造空間
をサンプルすることができる，拡張アンサンブル法の使用が必須であり，我々の新手法の開発によって，初めて現実 的な問題になったと言える。我々の結論は既存のどの力場も完璧なものはないというものであった。そこで，我々は 独自の力場パラメターの提案を我が国から初めておこなった（Sakae and Okamoto, Chem. Phys. Lett. (2003)）。我々 は更にタンパク質の折り畳みの遷移状態を詳しく調べることができる新しい拡張アンサンブル法を開発した（Berg, Noguchi and Okamoto, Phys. Rev. E (2003)）。拡張アンサンブル法では，平均値は求められるが，分子動力学の情報は 失われる。よって，我々は，モードカプリング理論と拡張アンサンブル法を組み合わせた，新しい分子動力学の研究 手法も提案した（La Penna, Mitsutake, Masuya and Okamoto, Chem. Phys. Lett. (2003)）。また，我々は独自に開発した拡 張アンサンブル法である，レプリカ交換アンブレラサンプリング法を D NA の塩基スタッキングの自由エネルギー 計算に適用し，実験と一致する結果を得た（Murata, Sugita and Okamoto, Chem. Phys. Lett. (2004)）。最後に，拡張アン サンブル法による膜タンパク質の立体構造予測法の提案も行った（Kokubo and Okamoto, Chem. Phys. Lett. (2004)）。 b) 生体分子の系以外にもエネルギー極小状態が多数存在する複雑系では，拡張アンサンブル法の適用が有効である。 本年度は，これまで我々が扱ってきた拡張アンサンブル法がカノニカルアンサンブル（定積定温アンサンブル）を元 にした手法であるのに対し，マルチカノニカル法を定圧定温アンサンブルに拡張し，唯一回のシミュレーションの 結果から，任意の圧力及び温度における定圧定温アンサンブル平均が得られる，新しい拡張アンサンブル法（マルチ バーリック・マルチサーマル法）を開発した（Okumura and Okamoto, Chem. Phys. Lett. (2004)）。この手法はタンパク 質の高圧変性の研究に応用できる。
B -1) 学術論文
A. MITSUTAKE, Y. SUGITA and Y. OKAMOTO, “Replica-Exchange Multicanonical and Multicanonical Replica-Exchange Monte Carlo Simulations of Peptides. I. Formulation and Benchmark Test,” J. Chem. Phys. 118, 6664–6675 (2003). A. MITSUTAKE, Y. SUGITA and Y. OKAMOTO, “Replica-Exchange Multicanonical and Multicanonical Replica-Exchange Monte Carlo Simulations of Peptides. II. Application to a More Complex System,” J. Chem. Phys. 118, 6676–6688 (2003). G. LA PENNA, A. MITSUTAKE, M. MASUYA and Y. OKAMOTO, “Molecular Dynamics of C-Peptide of Ribonuclease A Studied by Replica-Exchange Monte Carlo Method and Diffusion Theory,” Chem. Phys. Lett. 380, 609–619 (2003). B. A. BERG, H. NOGUCHI and Y. OKAMOTO, “Multi-Overlap Simulations for Transitions between Reference Configurations,” Phys. Rev. E 68, 036126 (11 pages) (2003).
Y. SAKAE and Y. OKAMOTO, “Optimization of Protein Force-Field Parameters with the Protein Data Bank,” Chem. Phys. Lett. 382, 626–636 (2003).
B -2) 国際会議のプロシーディングス
Y. OKAMOTO, “Metropolis algorithms in generalized ensemble,” in The Monte Carlo Method in the Physical Sciences: Celebrating the 50th Anniversary of the Metropolis Algorithm, J. E. Gubernatis, Ed., American Institute of Physics, Melville, pp. 248–260 (2003).
T. HIROYASU, M. MIKI, S. OGURA, K. AOI, T. YOSHIDA and Y. OKAMOTO, “Atom-level simulations of protein folding,” Proceedings of the 7th World Multiconference on Systemics, Cybernetics and Informatics (SCI2003), XIV, 117–122 (2003).
B -3) 総説、著書
岡本祐幸 , 「第一原理からのタンパク質立体構造予測の発展」, 特集「物理科学この１年」内 , パリティ 18, 54–56 (2003). 岡本祐幸, 「蛋白質折り畳みの計算機シミュレーション」, 特集「数学者のための分子生物学入門―新しい数学を造ろう」 内 , 物性研究 81, 93–119 (2003).
奥村久士、米沢富美子, 「液体水銀の金属−非金属転移領域における体積粘性率異常の理論」, 固体物理 38, 57–64 (2003).
B -4) 招待講演
岡本祐幸, 「タンパク質折りたたみの計算機シミュレーション」, 数学者のための分子生物学入門セミナー, 京阪奈, 2003年 1月 .
岡本祐幸, 「拡張アンサンブル法による蛋白質の立体構造予測」, 日本原子力研究所生体分子シミュレーションワークショッ プ , 京都 , 2003年 1 月 .
岡本祐幸 , 「タンパク質フォールディングの計算物理学：１次構造から３次構造へ」, 日本物理学会領域 12シンポジウム「生 物と物理学の架け橋―田中豊一が見たもの―」, 仙台 , 2003 年 3 月 .
奥村久士 , 「２体分布関数と粒子間ポテンシャル用いて書き表した体積粘性率の新しい表式」, 名工大物理工学セミナー, 名古屋 , 2003 年 4 月 .
Y. OKAMOTO, “Free energy landscape of protein folding studied by generalized-ensemble algorithms,” CECAM Workshop: Metastability and Landscapes in Complex Systems, Lyon (France), May 2003.
岡本祐幸 , 「拡張アンサンブルシミュレーションによる膜タンパク質の立体構造予測」, 分子研研究会「ロドプシンの分子科 学」, 岡崎 , 2003 年 5 月 .
Y. OKAMOTO, “Metropolis algorithms in generalized ensemble,” The Monte Carlo Method in the Physical Sciences: Celebrating the 50th Anniversary of the Metropolis Algorithm, Los Alamos, New Mexico (U. S. A. ), June 2003.
Y. OKAMOTO, “Protein folding dynamics studied by generalized-ensemble simulations,” Telluride Workshop: Protein Dynamics, Telluride, Colorado (U. S. A. ), July 2003.
Y. OKAMOTO, “Atom-level simulations of protein folding,” The 7th World Multiconference on Systemics, Cybernetics and Informatics (SCI2003) “Atom-Level Informatics, Simulations and Experiments in Materials Designing,” Orlando, Florida (U. S. A. ), July 2003.
奥村久士, 「体積粘性率の新しい表式 液体水銀の金属−非金属転移領域における異常の解釈」, 第１５回液体の化学夏 の学校 , 京都 , 2003年 8月 .
岡本祐幸, 「拡張アンサンブル法によるタンパク質の立体構造予測」, 電子情報技術産業協会バイオインフォマティックス専 門委員会 , 東京 , 2003年 9 月 .
岡本祐幸 , 「タンパク質折り畳みの統計力学」, 第８回久保記念シンポジウム, 東京 , 2003年 9 月 .
奥村久士, 「新しい拡張アンサンブルシミュレーション―マルチバーリックマルチサーマルモンテカルロ法―」, 東大工 学部 宮下・伊藤・藤堂研合同セミナー, 東京 , 2003年 10月 .
岡本祐幸 , 「タンパク質折り畳みの計算機シミュレーション」, 第７回シミュレーション・サイエンス・シンポジウム, 土岐 , 2003 年 10 月 .
岡本祐幸 , 「分子シミュレーションによる生体系の理論的研究」, 電力中央研究所セミナー , 我孫子 , 2003年 10月 .
岡本祐幸, 「拡張アンサンブル法によるタンパク質の折り畳みシミュレーション」, 早稲田大学理工学部セミナー, 東京, 2003 年 10月 .
B -7) 学会および社会的活動 学会誌編集委員
生物物理 会誌編集委員会委員 (2001-2002). 物性研究 各地編集委員 (2002- ).
Journal of Molecular Graphics and Modelling, International Editorial Board (1998-2000). Molecular Simulation, Editorial Board (1999- ).
科学研究費の研究代表者、班長等
日本学術振興会未来開拓学術研究推進事業「第一原理からのタンパク質の立体構造予測シミュレーション法の開発」, プロジェクトリーダー (1998-2003).
高校生対象の講演
「生体分子の計算機シミュレーション」, 愛知県立一宮高校スーパーサイエンスプログラム, 岡崎 , 2003年 11月 .
B -8) 他大学での講義、客員
SOKENDAI Lectures, “Protein folding simulations by generalized-ensemble algorithms,” Singapore National University, Singapore, January 11, 2003.
C ) 研究活動の課題と展望
我々が開発した新しい拡張アンサンブル法（特に，レプリカ交換マルチカノニカル法とマルチカノニカルレプリカ交換法）が 特に有効であることを示すことができた。そして，これらを駆使して，A MB E R やC HA R MMなどの生体高分子系における標 準的なエネルギー関数（力場）の是否の判定をすることができた。また，我々は，より精度の高いエネルギー関数を独自に開 発することにも成功したが，この研究は更に進める必要がある。特にベストと思われるものができ次第，ホームページで公開 し，世界の研究者に使ってもらう予定である。また，これまで水溶性タンパク質の立体構造予測に重点を置いてきたが，本年 度，膜タンパク質の立体構造予測に拡張アンサンブル法の適用が有効なことを示すことができた。特に，２００３年４月にヒトゲ ノム計画による塩基配列解読完了宣言がなされたが，これらのゲノム計画によって分かったことは，いろいろな生物体にお いて，遺伝子の約４分の１が膜タンパク質であるということである。現在，実験で決定されたタンパク質の立体構造が２万個以 上Protein D ata B ank（PD B ）に登録されているが，そのうち，膜タンパク質の立体構造は数十個に過ぎない。膜タンパク質は 結晶化が難しく実験によって，立体構造を決めるのが極めて困難であるからである。よって，拡張アンサンブルシミュレーショ ンによる膜タンパク質の立体構造予測はこれから重要性を増して行くであろう。
分子基礎理論第二研究部門
中 村 宏 樹（教授）
A -1)専門領域：化学物理理論、化学反応動力学論
A -2)研究課題： a) 化学反応の動力学
b)非断熱遷移の基礎理論の構築と応用 c) 化学動力学のレーザー制御
d)多次元トンネル理論の構築と応用 e) 超励起分子の特性と動力学 f) 分子機能の開発を目指して
A -3)研究活動の概略と主な成果
a) 化学反応の動力学：実際の化学反応において電子状態の変化する電子的に非断熱な化学反応が重要であることに鑑 み，かかる反応を取り扱う為の理論的手法の開発と具体的応用を進めている。正確な量子力学的な計算については， ３原子系でポテンシャルエネルギー曲面の交差領域で透熱表現を仮定する範囲内で（計算の規模が大きくなりはす るが），原理的には計算可能なコードを既に開発している。一方，大次元系にも適用出来る様にする為に古典および 半古典力学的な手法の開発を進めている｡古典的手法では，T S H（ T rajectory S urface Hopping）法を用い，ホップの確 率計算には我々独自の Z hu-Nakamura理論での確率表式を用いる。古典的に許されない遷移をも正しく取り扱える ので，例えば，L andau-Z ener公式を用いた場合よりも遥かに良いことが既に判っている｡L andau-Z ener型交差を持つ D H2
+
でそれを既に示し，現在は非断熱トンネル（Nonadiabatic T unnel）型遷移を有する系についての計算を実行して いる。以下に述べるトンネル理論をも取り入れて大きな系にも適用出来る手法を開発している｡更に，位相の効果を も含めた半古典力学的手法の開発も進めている｡これには，凍結波束の伝播を用いる Herman-K luk 理論に位相まで 含めたZ hu-Nakamura理論を組み込んでいる。１−２次元系でのテストを終えよく働くことを確認している｡量子力 学的計算では，O(1D )HC l 系の計算を進めている｡
b)非断熱遷移の基礎理論の構築と応用：Z hu-Nakamura理論の応用については，上記の反応を含めて以下で述べるc)，f) 等でも進めている｡基礎理論の開発では，先ず，空間座標と時間の両方に依存する（偏微分方程式の）問題に対する高 エネルギー近似での公式を求めることに成功した｡ポテンシャルが空間座標の一次関数で周波数が時間の一次関数 で変化するレーザー場が存在する場合の解である｡周波数が二次関数の場合にも解が得られることが判っている。 第二の課題は，レーザー場中での各種原子衝突過程の問題を高及び低周波数近似の下で取り扱うものである｡第三 の課題はボーズ・アインシュタイン凝縮系における原子の光会合過程（時間依存の非線形非断熱遷移問題になる）の 問題である｡非線形結合微分方程式の問題で，レーザー周波数を旨く制御することによって分子生成の効率を制御 する基礎理論が構築されうる｡
c) 化学動力学のレーザー制御：レーザーによる分子過程の制御に関する我々独自のアイディアは①フロケ表示，つま り光子の衣を着た状態の表現を用いること，及び，②この表現での非断熱遷移を制御することである。「着衣状態表
現での解析」はレーザーの周波数が遠赤外の様な低周波数の場合を除いて大変良い近似で成立する。また，状態遷移 は非断熱遷移を介して最も有効に起こる訳で，それに対する解析的基本理論を我々は有しており，しかも，周期チャ
―プによる完全制御法の理論を構築している。これらの理論は，当然ながら，エネルギー準位に対するものなので， 波束の運動を扱う時には，これら解析理論に基づいて得られる制御パラメーターを第一近似として更に改良を加え る努力を行う｡この基本姿勢に基づき，波束の運動，即ち，断熱ポテンシャル上の動きと非断熱遷移によるポテンシャ ルエネルギー曲面間の乗り移りを制御することを考える｡先ず，２原子分子でそのポテンシャルエネルギー曲線の 形や質量の違いによって色々な制御方法を提唱しようとしている。さらに，多次元系への拡張を可能にする為に最 適制御理論を古典軌道に基礎を置く半古典力学で再構成し直すことを考えている｡量子力学的な最適制御理論では， 数値計算がどうしてもグリッド法に基づくこととなり，実際の計算が２次元系までに限られている｡我々は，半古典 理論を用いてこの壁を破りたいと思っている。既に基本的定式化は終えている｡
d) 多次元トンネル理論の構築と応用：トンネルには，①二重井戸におけるエネルギー分裂，②トンネルを介しての崩壊，
③反応におけるトンネル，の三種類の現象がある。いずれについても，３次元以上の多次元系に適用出来る正確な理 論はなかったが，我々は，最初の二つの問題に対してインスタントン理論を実用性のあるものに定式化し直す事に よって高い精度で高次元系を扱い得る理論の開発に成功した。インスタントンと称する周期軌道を効率良く見出す 手法の開発と座標変換に対して不変な定式が基本である｡９原子系のマロンアルデヒド分子の水素結合を介しての 水素移動に適用し実験との極めて良い一致を得ている｡高次元系で高精度の量子化学計算を限られた空間配置で行 うだけで精度良くエネルギー分裂等を評価する手法をも開発した｡振動励起状態に対する理論も構築した｡一方，古 典あるいは半古典論的な反応動力学にとって重要な問題の一つはトンネル効果を如何に取り入れるかである。それ には先ず，古典軌道上で火線（１次元系での転回点に相当）を効率良く見出す方法の開発が必要である。我々は，軌道 がカオス的に振舞う場合であっても火線を効率良く見出す手法を開発した｡この火線からインスタントン的トンネ ル軌道を走らせることによってトンネル効果を反応に取り入れることが出来る。
e) 超励起分子の動力学：二種類の超励起状態（多電子或いは内殻電子が励起した状態（第一種），及び振動回転励起した リュードベリー状態（第二種））が関与する動力学を律する重要なパラメータとして，夫々のポテンシャルエネルギー 曲線以外に電子結合（第一種の状態と連続状態の結合）と量子欠損がある。我々は最近，２原子分子の核間距離に依 存した量子欠損を電荷分布から解析的に評価する理論を構築することに成功した｡核間距離依存性を実験から求め ることは容易ではないので，有効な理論となり得る。
f) 分子機能の開発を目指して：我々は非断熱トンネル型の遷移で完全反射現象が起こることを見出しているが，これ を積極的に利用することによって分子に機能を持たせることが出来るはずであると考えている。また，フォトクロ ミズム等の機能においても非断熱遷移が重要な役割をしている筈であり，非断熱遷移を積極的に制御することに よって，分子の新しい機能を開発できる筈であると考えている｡具体的分子を例にとって電子状態の量子化学計算 をも行いつつ各種可能性を追求している｡
B -1) 学術論文
A. KONDORSKIY and H. NAKAMURA, “Photodissociation of H2+ and HD+ in an Intense Laser Field,” Phys. Rev. A 66, 053412 (8 pages) (2002).
G. V. MIL’NIKOV and H. NAKAMURA, “Calculation of Resonances in the dtµ Molecule by the R-Matrix Method,” Phys.Rev. A 67, 034501 (4 pages) (2003).
G. V. MIL’NIKOV, K. YAGI, T. TAKETSUGU, H. NAKAMURA and K. HIRAO, “Tunneling Splitting in Polyatomic Molecules: Application to Malonaldehyde,” J. Chem. Phys. 119, 10–13 (2003).
H. NAKAMURA and E. A. SOLOVJEV, “Semiclassical Theory for the Quantum Defect Function of Diatomic Molecules,” J. Phys. B 36, 3697–3705 (2003).
H. NAKAMURA, “Analytic Solution to Wave Packet Dynamics in a Laser Field: the Case of Linear Chirp,” Chem. Phys. 295, 269–273 (2003).
B -3) 総説、著書
L. PICHL, M. HIYAMA and H. NAKAMURA, “Analytical Treatment of the K-Matrix Integral Equation in the Dynamics of Superexcited Molecules,” in Dissociative Recombination of Molecular Ions with Electrons, S. L. Guberman, Ed., Kluwer Academic/Plenum Publ., p.167–176 (2003).
H. NAKAMURA, “Nonadiabatic Transitions and Chemical Dynamics,” in Current Developments in Atomic, Molecular, and Chemical Physics with Applications, Man Mohan, Ed., Kluwer Academic/Plenum Publishers; New York, p.71–77 (2002).
B -4) 招待講演
中村宏樹 , 「レーザーで反応を制御する」, 日本学術会議講演会「原子分子光科学の最前線」, 東大理化学 , 2003年 6 月 . 中村宏樹 , 「レーザーによる反応制御の可能性を操る」, 特定領域「強光子場分子制御」領域セミナーにて招待講演 , 東大 理化 , 2003 年 5 月 .
H. NAKAMURA, “Quantum and Semiclassical Dynamics of Chemical Reactions, Multi-Dimensional Tunneling, and Laser Control of Molecular Processes,” Workshop on “Quantum Reactive Scattering,” San Lorenzo de El Escorial (Spain), June 2003.
H. NAKAMURA, “Theory of Tunneling Splitting in Polyatomic Molecules: Full Dimensional Treatment of Hydrogen Transfer in Malonaldehyde,” Multidimensional Quantum Reaction Dynamics 2003, Berlin, July 2003.
H. NAKAMURA, “New Developments in the Theory of Nonadiabatic Transitions: Zhu-Nakamura Theory and Some Others,” Moscow–Chernogolovka Workshop on Nonadiabatic Transitions in Quantum Mechamics, Moscow, Chernogolovka, August 2003.
H. NAKAMURA, “Applications of the Zhu-Nakamura Theory to Electronically Nonadiabatic Chemical Dynamics,” G. Billing Memorial Symp., Copenhagen, August 2003.
H. NAKAMURA, “Semiclassical Theory of Laser Control of Molecular Processes: Importance of Nonadiabatic Transition,” Intern. Symp. on Ultrafast Intense Laser Science 2: Propagation and Interaction, Lac Delage, Quebec (Canada), September 2003.
B -6) 受賞、表彰
中村宏樹 , 中日文化賞 (2000).
B -7) 学会および社会的活動 学協会役員、委員
原子衝突研究協会委員 (1981-1994). 学会の組織委員
IC PE A C（原子衝突物理国際会議）第９回組織委員会 , 経理担当 (1979). IC PE A C（第１７回及び第１８回）全体会議委員 (1991, 1993).
IC PE A C（第２１回）準備委員会委員 , 運営委員会委員 (1999).
A IS A MP（アジア原子分子物理国際シンポジウム）A dvisory committeeメンバー (1997, 2002). Pacifichem 2000 シンポジウム組織者 (2000).
Workshop on Nonadiabatic Transitions in Quantum Mechanics. Internat. Advisory Committee Member (Moscow- Chernogolovka, August 2003).
文部科学省、学術振興会等の役員等
学術審議会専門委員 (1991-1995, 1998-2002, 2002- ). 理化学研究所基礎科学特別研究員審査委員 (2003- ). 学会誌編集委員
Computer Physics Communication, Specialist editor (1986- ).
Journal of Theoretical and Computational Chemistry, Executive editor (2001- ).
科学研究費の研究代表者、班長等
特定領域研究計画班代表者 (1999-2001). 基盤研究代表者 (1998-2000, 2001-2003). 特別推進研究代表者 (2003- ).
岡崎高校スーパーサイエンスハイスクール活動支援 (2002- ) 分子研総括責任者 .
講演「学問創造への挑戦―未来をになう皆さんへ」.
B -8) 他大学での講議、客員
ウォータールー大学応用数学科 , ウォータールー（カナダ）, -2003年 6月 .
分子基礎理論第四研究部門
平 田 文 男（教授）
A -1)専門領域：理論化学、溶液化学
A -2)研究課題：
a) 溶液内分子の電子状態に対する溶媒効果と化学反応の理論 b)溶液中の集団的密度揺らぎと非平衡化学過程
c) 生体高分子の溶媒和構造の安定性に関する研究 d)界面における液体の統計力学
A -3)研究活動の概略と主な成果
a) 溶液内分子の電子状態に対する溶媒効果と化学反応の理論：溶液中に存在する分子の電子状態は溶媒からの反作用 場を受けて気相中とは大きく異なり，従って，分子の反応性も違ってくる。われわれは以前にこの反作用場を液体の 積分方程式理論によって決定する方法（R ISM-SC F 法）を提案している。この理論を使って２００３年度に行った研究の 主な成果を以下にまとめる。
(i)溶液内有機化合物のイオン過程に伴う電子構造と溶媒の再配置：溶液内電子移動反応においてはいわゆる溶媒の再配 置（向）エネルギー（あるいは非平衡自由エネルギー）が本質的役割を演じることはよく知られている。われわれは以前にい わゆるマーカスの自由エネルギー曲面を分子レベルで記述する方法をR IS M理論と熱力学サイクル（R IS M-T C 法）を使っ て提案している。[J. Phys. Chem. 99, 10526 (1995)]一方，K oboriらは電子スピン共鳴の結合定数（J）を解析することによっ て溶媒の再配置エネルギーを決定する実験的方法を提案し，N,N-deimethylaniline（D MA ）を始めとするいくつかの溶質の 電子移動反応に伴う溶媒の再配置エネルギーを様々な溶媒中で決定している[J. Am. Chem. Soc. 123, 9722 (2001)]。本研 究ではR ISM-SC F 理論とR ISM-T C 法を組み合わせた新しい方法によりアセトニトリル中のD MA のイオン化過程に伴う溶媒 の再配置エネルギーおよび電子構造の変化を求め，それらの間の相関について考察を行った。[J. Chem. Phys. 119, 2753 (2003)に既報]
(ii)溶質−溶媒間の相互作用を量子論的に記述した溶液内電子状態理論の提案：溶液中における種々の化学過程におい て，電子に起因する現象は枚挙に暇がない。例えば，溶液中の化学反応はその好例であり，科学の根幹をなす非常に重要 な現象である。近年，このような現象を理論的に取り扱うために，連続誘電体モデル，QM/MM 法，R IS M-S C F 法といった溶 液中分子の電子状態理論が提案されているが，既存の方法の多くは，分子間相互作用に古典的相互作用を用いるという 理論的な弱点を持っている。そこで我々は，気相中の電子状態理論と同等の意味での非経験的な溶液中分子の電子状態 理論の構築の第一歩として，単純液体を対象に新規な電子状態理論を提案した [T. Yamazaki, H. Sato and F. Hirata, J. Chem. Phys. 119, 6663 (2003)]。本理論では，上述の理論とは異なり，溶質−溶媒間相互作用は量子論的に記述される。今 回，その理論を拡張し，水分子等の分子性溶媒を扱うことのできる理論を構築した。量子論的な近距離力に基づく反発が強 く出過ぎるという問題はあるものの，水溶液中の単原子分子の電子状態に関して一定の成果を挙げており，現在投稿準備
中である。
b)溶液中の集団的密度揺らぎと非平衡化学過程：われわれは昨年までの研究において，液体の非平衡過程を記述する
上で相互作用点モデルが有効であることを示し，そのモデルによって液体中の集団的密度揺らぎ（集団励起）を取り 出す方法を提案してきた。さらに，その理論に基づき溶液内の化学種のダイナミックス（位置の移動，電子状態，構造 変化）をそれらの変化に対する溶媒の集団的密度揺らぎの応答として記述する理論を展開しつつある。この分野の 研究の主な成果は以下のとおりである。
(i)溶媒和ダイナミクスの溶質構造依存性：S S S V 理論による研究：われわれは以前にR IS M 理論と一般化ランジェヴァン方 程式を結合し，過減衰領域で溶媒のダイナミクスを記述する理論，Site-site Smoluchowski-V lasov（SSSV ）理論を提案し，水 の構造ダイナミクス[J. Chem. Phys. 96, 4619 (1992)]，イオンの荷電状態の変化に誘起された溶媒のダイナミクス[J. Mol. Liq. 65/66, 15 (1995)]，アセトニトリル中ベンゾニトリルの垂直励起状態から平衡励起状態への緩和に伴う溶媒および電子構造 の緩和ダイナミクス[J. Chem. Phys. 110, 11423 (1999)]，などの問題を取り扱った。本研究では溶質の幾何形状および電荷 分布が異なる様々なモデル系にSSSV 理論を適用し，溶質の幾何形状や電荷分布が溶媒和ダイナミクスに及ぼす影響を考 察した。[J. Chem. Phys. 118, 2279 (2003)に既報]
(ii)水中の溶質の回転緩和に対する圧力効果：溶液内の一個の溶質分子のダイナミクスは溶液内化学反応ダイナミクスを 定式化する上で重要なステップとなる。溶液内化学反応ダイナミクスの問題ではこれまで伝統的に溶媒を連続媒体で塗る つぶす乱暴なモデルが使われている。例えば，電子移動反応のマーカス理論やいわゆるクラマース理論がその典型例で ある。このような理論では溶質は量子力学を使ってまさに原子レベルで表現するのに対して，溶媒はマクロな誘電率や粘性 で特徴づけるため，その記述は空間的にも時間的にもその分解能に著しい違いがある。一方，現在の実験は溶媒の構造や ダイナミクスを原子レベルの解像度で測定することができる段階に到達しており，実験的分解能に対応可能な理論の構築 はまさに焦眉の課題である。本研究ではそのような定式化への第一ステップとして，R IS M理論とモード結合理論を組み合 わせた方法により，溶液内の一個の溶質の回転ダイナミクスを記述する理論を定式化し，その理論により水中におけるアセ トニトリルおよびメタノールの回転緩和速度に対する圧力効果を解析した。[J. Mol. Liq.に投稿中]
c) 生体高分子の溶媒和構造の安定性に関する研究：本研究課題の最終目的は第一原理すなわち分子間相互作用に関す る情報のみから出発して蛋白質の立体構造を予測することである。蛋白質の立体構造予測（すなわちフォールデイ ング）には二つの要素がある。そのひとつは広い構造空間をサンプルするための効果的なアルゴリズムであり，他は 蛋白質の構造安定性を評価する問題である。蛋白質の安定性はそれが置かれている環境すなわち熱力学的条件に よって完全に規定される。この熱力学的条件には溶媒の化学組成（溶媒の種類および共存溶質の濃度），温度，圧力な どが含まれる。本研究では蛋白質の構造安定性に対して熱力学的条件が与える影響を分子レベルで明らかにする目 的で，その素過程として，アミノ酸やペプチドおよび疎水分子の水和現象を分子性液体の統計力学（R ISM理論）に基 づき解析している。これらの解析は蛋白質の安定性に関わる物理的要因を分子レベルで解明するだけでなく，今後， 蛋白質のフォールデイングを実際に実行するうえで重要となる溶媒和自由エネルギーを計算するための方法論的 基礎を与えるものである。
(i)分子内自由度をもつ溶質の部分モル容積および部分モル圧縮率の理論：これまで我々のグループではR IS M理論に基 づき，溶液の部分モル容積（PMV ）や部分モル圧縮率（PMC ）を求める理論を開発し，水溶液中の20個のアミノ酸の部分モ ル容積を，ほぼ，完璧に再現すると同時に，ペプチドのヘリックス−コイル転移に基づく部分モル容積変化を定性的に求め ること成功している。しかしながら，これらの計算では溶質分子の幾何構造が固定していることを仮定しており，溶質の構造 揺らぎの影響は考慮に入れていない。実際には自由度の大きな分子の場合，その構造自身も熱運動によってある程度揺ら いでいることが予想されるため，この揺らぎが部分モル容積や部分モル圧縮率にどのような影響を及ぼすかという問題はこ れらの実験結果から溶媒和に関する情報を得る上で本質的な問題のひとつである。本研究では分子内自由度をもつ溶質
の部分モル容積および部分モル圧縮率を取り扱う一般的な統計力学理論を提案すると同時に，これをR ISM理論と結合し て，水中のブタンの異性化反応に関するモデル計算を行った。その結果，下記の諸点を明らかにした。①平衡系のPMV は 個々の異性体のPMV の統計平均で表現される。②平衡系のPMC は個々の異性体のPMC の統計平均と構造変化（反応） 座標に沿った PMC の平均自乗揺らぎの和で表現できる。[J. Chem. Phys. 119, 5623 (2003)に既報]
d)界面における液体の統計力学：
(i)R eplica-R IS M 理論による電気二重層キャパシタの研究：電気二重層キャパシタは，多孔質炭素電極の細孔に形成され る電気二重層に電荷を蓄えるデバイスであり，通常のコンデンサーの100万倍から1000万倍の電気容量を示す。安価で無 害なことから，多孔質炭素電極を使った高容量な蓄電デバイスの発展が期待されている。最近，通常の活性炭（A C ）と炭化 したポリ塩化ヴィニルデン（PV D C ）の電気二重層容量についての比較実験が行われ，PV D Cの方がより高容量であること が示された。そこで，本研究では，以前に，当研究室で開発したR eplica-R IS M理論を用いて，電極の細孔構造と電気二重 層容量の関係について調べ，二重層容量の決定因子について考察した。[Chem. Phys. Lett. 378, 638 (2003)に既報] (ii)気液界面および気液相転移の統計力学理論：気体と液体あるいは液体と液体など流体間の界面は平均の密度（あるい は濃度）が空間的に変化することから，一様な流体を対象とする通常の液体統計力学を適用することができない。このような 問題ではいわゆる二体相関関数が一体の相関関数（平均密度）に依存するだけでなく，逆に，一体相関関数が二体相関関 数に依存するからである。したがって，流体界面の構造を明らかにするためにはにこれら二つの相関関数に関する連立積 分方程式を自己無撞着的に解法することが必要となる。このような理論の基本的枠組みはすでに密度汎関数理論として提 案されているが，この理論を実際に気液界面や液液界面の問題に適用した例はない。われわれは本研究において，L ovett， Mou，B uff およびW ertheimによって提案された密度汎関数理論（L MB W 理論）を 気液界面および気液相転移の問題に適 用し，気液共存曲線を求めると同時に，界面における液体構造を明らかにした。この研究における最も重要な結論は気液界 面における二体相関関数が臨界点におけるそれに酷似した非常に長い相関長を示すことであり，このことは気液界面にお いてマクロなスケールの密度揺らぎが存在することを意味している。これまで，実験的にはそのような密度揺らぎが知られて おり，「リプロン」とか「キャピラリウエーブ」という言葉で表現されていたが，このようなマクロな密度揺らぎを第一原理的に再 現したのはこの理論が最初である。[J. Chem. Phys.に投稿中]
(iii)ブタノール−水系の圧縮率：アルコール−水混合系は人類にとってもっとも身近な溶液である。しかしながら，その構造や 物性に関する理解は２１世紀に至っても未だ「渾沌」と言っても過言ではない。例えば，アルコールと３級ブタノールはあらゆる 濃度で混合し，見掛上，水相とアルコール相への相分離は観測できない。しかしながら，ある温度，圧力では相分離にも近い 大きな濃度揺らぎを示すことが知られている。このような事実から，多くの研究者が水−アルコール中にはミクロスケールでは あるがかなり大きな濃度の揺らぎが存在すると考え，これを「ミクロクラスター」とか「疑似相分離」という言葉で表現してきた。 最近では水−プロパノール系にイオンを添加した場合に水相とアルコール相への相分離が観測された。（新潟大，三沢）こ れは水−アルコール系が熱力学的な「準安定」状態に近いような混合相を作っていることを示唆している。それでは，その「ミ クロクラスター」や「疑似相分離」なるものの構造的実態は一体何か？ われわれは，最近，R ISM理論に基づき水−エタノー ル系および水−ブタノール系の溶液構造に関する解析を行った。さらに，水−ブタノール系に関してはその熱力学的安定性
（圧縮率および濃度揺らぎ）を解析した結果，下記のような結論を得ている。
①純水ではこの温度領域において圧縮率が温度上昇とともに減少する。この挙動は水中のバルキーな（圧縮され易い）氷 様構造が温度とともに減少するという描像で大雑把に説明できる。②純ブタノール系では圧縮率が温度とともに単調に増加 するが，これは通常の液体に共通の挙動であり，熱膨張によって圧縮できる空間が増大するという物理で説明することがで きる。③圧縮率の濃度依存性はこの温度領域において極小値を示す。この極小値の出現は圧縮率の濃度依存性に二つの
相反する傾向が存在し，それらがある濃度で拮抗することを示している。それらの傾向のうちのひとつはブタノールの高濃 度領域における濃度増加に伴う圧縮率の増大であり，これは理想混合モデルで説明できる。すなわち，混合系の圧縮率を
（χT = (x–1)χwater + xχButanol）のように水の圧縮率とブタノールのそれの成分比の重みをつけた単純な重ね合わせとみ
なせば，水に比べて大きな圧縮率をもつブタノールの寄与が濃度の増加とともに支配的となり，高濃度領域では系全体の圧 縮率は濃度とともに増加する。一方，低濃度領域で濃度とともに圧縮率が減少する傾向は，そのような単純な描像では説明 できない。この傾向を説明する上で，水−ブタノール系の液体構造に関する上述の結果が有効である。すなわち，上述のよ うに，低濃度領域においてブタノールはその水酸基を介して水の水素結合ネットワークの中に組み込まれ，そのアルキル基 は水の氷様構造の隙間の中に入り込んでいる。すなわち，アルキル基が水の圧縮可能な空間を充たすことにより圧縮率を さげているのである。[J. Theor. Comput. Chem. 2, 193–203 (2003)に既報]
B -1) 学術論文
T. SUMI and F. HIRATA, “A Density-Functional Theory for Polymer Liquids Based on Interaction Site Model,” J. Chem. Phys. 118, 2431–2442 (2003).
K. NISHIYAMA, F. HIRATA and T. OKADA, “Solute-Structure Dependence of Solvation Dynamics Studied by Reference Interaction-Site Model Theory,” J. Chem. Phys. 118, 2279 (2003).
T. YAMAGUCHI, S. -H. CHONG and F. HIRATA, “Dielectric Relaxation Spectrum of Water Studied by the Site-Site Generalized Langevin/modified Mode-Coupling Theory,” Mol. Phys. 101, 1211–1220 (2003).
T. YAMAGUCHI, S. -H. CHONG and F. HIRATA, “Theoretical Study on the Molecular Motion of Liquid Water under High Pressure,” J. Chem. Phys. 119, 1021 (2003).
T. IMAI and F. HIRATA, “Partial Molar Volume and Compressibility of a Molecule with Internal Degrees of Freedom,” J. Chem. Phys. 119, 5623–5631 (2003).
H. SATO, Y. KOBORI, S. TERO-KUBOTA and F. HIRATA, “Theoretical Study on Electronic and Solvent Reorganization Associated with a Charging Process of Organic Compounds: I. Molecular and Atomic Level Description of Solvent Reorganization,” J. Chem. Phys. 119, 2753–2760 (2003).
I. OMELYAN, A. KOVALENKO and F. HIRATA, “Compressibility of tert-Butyl Alcohol–Water Mixture: the RISM Theory,” J. Theor. Comput. Chem. 2, 193–203 (2003).
A. SETHIA, E. R. BITTNER and F. HIRATA, “Interplay between the Repulsive and Attractive Interaction and the Spacial Dimesionality of an Excess Electron in a Simple Fluid,” J. Theor. Comput. Chem. 2, 129–138 (2003).
A. TANIMURA, A. KOVALENKO and F. HIRATA, “Molecular Theory of an Electrochemical Double Layer in a Nanoporous Carbon Supercapacitor,” Chem. Phys. Lett. 378, 638–646 (2003).
T. YAMAZAKI, H. SATO and F. HIRATA, “A Quantum Solute-Solvent Interaction Using Spectral Representation Technique Applied to the Electronic Structure Theory in Solution,” J. Chem. Phys. 119, 6663–6670 (2003).
B -3) 総説、著書
F. HIRATA, “Molecular Theory of Solvation,” Kluwer-Springer Academic (2003).
